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Abstract 



o 
o 

CN ■ One of the promising sources of gravitational radiation is a binary system composed of compact stars. It is 

f-H , an important question how the rotation of the bodies and the eccentricity of the orbit affect the detectable 

■ signal. Here we present a method to evaluate the gravitational wave polarization states for inspiralling 
compact binaries with comparable mass. We consider eccentric orbits and the spin-orbit contribution in the 

, case of one spinning object up to 1.5 post-Newtonian order. For circular orbits our results are in agreement 

' with existing calculations. 

1 Introduction 
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During the inspiral and merger of compact stars intense gravitational radiation is produced, whose detection 
is expected by the current generation of gravitational wave observatories and more sensitive instruments, such 
as advanced LIGO^j and the LISA mission The different characteristics of the binary, e.g. rotation of the 
bodies and eccentricity of the orbit contribute to the emitted radiation. Having accurate information about 
the waveforms will significantly increase the possibility of detection and the precision with which the source 
QQ I parameters can be identified. 

. The major approximation tool for the description of the dynamics and the generated waveforms of a binary 

CO ' system is the post-Newtonian (PN) expansion f3J . Neglecting ultrarelativistic and extreme gravitational effects 
l/^ . in this weak field approximation the velocities and the gravitational potential involved are small and there is no 
' restriction on the mass ratio of the components. This approach precisely describes the motion of an inspiralling 
. binary up to the last stable circular orbit [H [S] . To explore the evolution of the binary in the merger or ringdown 
phases one has to use different methods [B] or numerical simulations to solve the Einstein equations [7] . 

The polarization states of the emitted gravitational waves was computed for quasi-circular [51 [SI [10] and 
' elliptic orbits [TTJ [T^l [13] • For spinning binaries the evaluation of the wave pattern has been done by several 
i authors [14[ [151 E] • In niost cases the waveform is expressed formally in terms of the dynamical quantities of 
' the motion. Here we give the expressions of the wave polarization states h+ and up to 1.5 relative PN order 
for eccentric orbits. Extending our former results |17| we discuss the effects of rotation for a binary system in 
which one of the bodies is spinning with spin vector S and the components have comparable masses mi and 

7712. 

In Ref. '18] the equations of motion for a comparable mass binary system are derived by the use of the 
Lagrangian formalism and an appropriate radial parameterization of the orbit was given |19j . We use these 
results to determine the features of the motion up to 1.5 PN order. In order to avoid long expressions we 
introduce the invariant and comoving coordinate systems which are fixed to the total and Newtonian angular 
momentum vectors, respectively. In Sec. 2 we evaluate the quantities which are necessary to obtain the general 
formulae of the detectable gravitational wave signals. In the first subsection we introduce the invariant and 
the comoving coordinate systems and give the formal expressions of the (N, p, q) triad which determines the 
relative orientation of the source and the observer. In the second subsection we describe the precession of the 
spin. In the third one, with the use of the constants of the motion, we give the equations for the radial and 
angular variables of the orbit linearly in spin. In Sec. 3 we recall that the detectable signal separates into 
and hx and calculate the formal expressions of the polarization states from the transverse-traceless tensor hl^j, 
which represents metric perturbations. Having a great importance, in Sec. 4 and 5 we investigate the circular 
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orbit case and the extreme mass ratio limit. Our results for circular orbits are in agreement with those of Kidder 
[15j . Sec. 6 and the Appendix contain our conclusions and the general expressions for the transverse-traceless 
tensor in our case. 

We use units in which c = G = 1 . 



2 Description of the motion 

2.1 Basic vectors in the invariant and comoving coordinate systems 

To describe the time dependence of the orbital elements we introduce a coordinate system which does not change 
in time. The z axis of this invariant system is fixed to the direction of the total angular momentum vector J. 
This vector is constant up to 2 PN order [TB] . We choose the x and y axes in a way that the vector representing 
the direction of the line of sight has the form N (sin 7, 0, cos 7) in this system, where 7 is the constant angle 
between J and N. 

As in the case of a Lense-Thirring system |17j we introduce the comoving coordinate system, in which the 
X and z axes are fixed to the separation vector r and the Newtonian angular momentum vector L^r = /ir x v, 
respectively. Here v is the relative velocity vector, ^ = 77111712 / M"^ is the reduced mass and M — mi + m2 is 
the total mass of the system. 

The transformation between the coordinate systems is described by Euler-angles |18j . A general vector u of 
the comoving system becomes 

u' = i?,(<I>)i?,(Oi?.(*)u (1) 

in the invariant one. Here l is the angle between J and Lat, <f> describes the precession of L^v over J and 4" 
represents the direction of the separation vector on the orbital plane, the plane perpendicular to Lat. This way 
the components of the separation vector in the invariant system are 

cos $ cos Vl/ — cos L sin <i> sin 
r I sin cos \1/ + cos l cos $ sin ^I* | . (2) 
sin i sin '5 

With the use of the Euler-angles the relative velocity vector, which is perpendicular to the Newtonian angular 
momentum, is expressed as 

f r \ 

(3) 




in the comoving coordinate system. Furthermore we decompose v as 

V = t;||n + v±m , (4) 

where n = r/r and m is the unit vector parallel to the y axis of the comoving system. 

To collect all the terms describing the spin effects of the rotating body we decompose the dynamical quantities 
into zeroth-order and linear terms in spin: 

r = rN+rs, W|| = ""11^ + , v±=vj_n + v±s 

* = + ^-S , $ = $Ar + $s , L = LN + LS . (5) 

For the calculation of the polarization states we have to determine the components of the orthonormal triad 
(N, p, q), where N is the direction of the line of sight and p is a vector perpendicular to N and Ljy. 

Although the evolution of the angles will be discussed later, we introduce an important result in advance, 
namely ln = iS]. It is convenient to introduce the angle T = + $ which is decomposed into zeroth order 
and linear terms in spin, cf. Eq. ([5]). In this case N has the form 

sin7C0sTAr — (sin7T5sinTjv — cos7sin$jv'-s) 
N= I — sin7sinTAr — (sin7T5CosTAr — cos7Cos$jvi5) | (6) 
cos 7 + sin7sin$jv''S 
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in the comoving system. Since p is a unit vector perpendicular to N and Lat it beco 



imes 



sinTjv + (T5 — cot7COS$Art5) cosTjv 
p= cosTat — (Ts — cot7Cos$Arts)sinTjv 





(7) 



in the comoving system and q = N x p. 

2.2 Spin-precession 

To describe the dynamics of the binary system first we have to determine the evolution of the spin vector S. In 
the invariant system let the angles a and (3 denote the direction of the spin, S = S'(sinQ! cos /?, sina sin/3, cos a). 
The dynamics of S is governed by the spin precession equations |20j 



S=(4 + 3C)^LArxS, 



(8) 



where C — "012/ mi. Up to 2PN order the magnitude S of the spin is constant. Since we are interested in the 
leading order spin effects we may replace Lat ^ J by inserting higher order terms in S. After substituting the 
components of the spin vector Eq. ^ leads to the following equations for a and (3: 



d = 



/3=(4 + 3C), , 



_J 

2r 



N 



The relative PN order of these angles can be determined with the use of T^r as zeroth order reference 



L/{rl) _ M 



~ e 



(9) 



(10) 



where e is the post-Newtonian parameter. We assume that the integration of the equations of motion does not 
change the order of the different quantities. Hence (3 can be decomposed as /3 = (3n + Ppn with /3jv/Tjv ~ 1 
and Ppn/^n e, and 



f3- 



N 



0, /3pjv = (4 + 3C)-^ 



(11) 



2.3 The equations of motion 

The length and the first component of the relative velocity vector are [18] 

2CLS 2 2£; 2M L2 2CSLS 2(2 + C)LS 



2E 



2M 
r 



2E 



2M 

r 



fir'' 



(12) 



respectively, where E is the energy, L denotes the length of the angular momentum vector L = L 
cosK = LS/LS. These quantities are constants of the motion yLSjj. Since v"^ = - 
of the relative velocity is 



N 



-'SO 



Vl_ = 



L L frs C-ELS 2/iLS\ 



Using Eq. ([3]) we get the equations of the motion for the angle T 



and 

the second component 
(13) 



Tc = — 



L 



/2rs (ELS 2^LS\ 
The remaining equations are determined by the total angular momentum J = 



L 



-'SO 



S, where 



-'SO 



77 



;2 + C)^[rx(rxS)]-|[vx(vxS)] 



(14) 



(15) 
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and rj = /i/M [151 HH]- After substituting the components of r, v and S and using the condition that J is a 
constant vector we get the following equations 



tssin<i>Ar = — (Acos^ + _B sin^ + Ccos/3jv — /3pjv sin/3jv) 



(Bcos^ — Asin^ — C sin/^Ar — /3pjv cos/3jv) 



tgCOS^AT = 

where ^ = ST^r — I3n and 

"(2 + C)/i vC{vIn-'^In) 



S sin a 
L 

S sin a 



A = 



2rN 



77C< (2 + C)m' 



2r 



N 



(16) 



(17) 



Although we cannot give a full description of the angular evolution we have determined all the quantities, 
namely the equations for the angle T and the products ls cos $jv and ls sin^pf we need to describe the evolution 
of the polarization states. 



3 The polarization states 

The signal h{t) of a laser-interferometric gravitational wave detector is decomposed into the polarization states 
h+{t) and h^it) [16], 



h{t) ^ F+h+{t) + F^h^it) 



(18) 



where and Fx are the so-called beam-pattern functions. The independent polarization states h^{t) and 
hx(t) are projected from the transverse-traceless tensor hlfj, representing metric perturbations as 



h+ = ^{ptPj - gj9j)/iTT > hx = ^{piQj + qiPjWj.j. 



(19) 



In the post-Newtonian approximation /i^j, can be decomposed as 



2^ 



TT 



(20) 



where D is the distance between the source and the observer. Q*^ denotes the quadrupole (or Newtonian) 
term, P^-^Q^^ , PQ^^ and P^-^Q^^ are corrections corresponding to higher PN orders, PQgp ^^"^ -^^^Q^so 
the spin-orbit terms [THl HI] . Since we are interested in the effects of rotation we keep the contributions linear 
in spin and the quadrupole term. The decomposition of the relative velocity vector, Eq. ([4]), gives a natural 
structure to hl^j, and the components can be described in a simple way, see Appendix A. 

To avoid complicated expressions the components of N, p, q, v and S are inserted formally. We decompose 
the relevant contributions to the polarization states /i+ and /ix as 



2^ 



(21) 



where 



l,5SO 



- —] (Pt- it) + '2-Vl^f{p^Py - q^Qy) + V'l{pt 



1^ [(qS)p, + (pS)(7.] 



\{iv^{i~(:)sM 



' 'Z^O + ^^^fSx (p,p-feq)]. 



(4 - bC)v^ 



[S X (pj;P - qx(\}\y - v^C[S X (pyp - qyC{% 



CS 



+ v±Ny (pa;q -I- q^p) + v^N^ipyq + qyp) 



(22) 
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and similarly 



l,5SO 



M 
r 



Pxqx + v±r{Pxqy + qxPy) + v\pyqy 



^ [(qS)fe - mpx] , 

2 f {2 + Or 
— < 6v±{l - OSzPxqx + 7; [S X {px(l + qxP)]a 



(4 - 5C)v^ 



[S X (p^q + qxP)]y - vj_C[S X (pyq + qyp)]j 



-N^ + vj_Ny ) (qj;q - p^p) + v^N^{qyC[ - pyp) 



(23) 



Eqs. (|22j|23p form the basics of our main results. With the use of these expressions one can evaluate the 
spin contributions to the polarization states. After parametrizing the orbit [21U19j the equations of the motion 
for the angles T, l, $0 and f3 can be integrated. With the substitution of these angles into Eqs. (I6][7]) the 
explicit parameter dependence of the (N,p,q) triad is determined. Using Eqs. (j22ll23p and neglecting quadratic 
or higher order spin terms one can investigate the effects of rotation on the detectable gravitational waveform. 

4 The circular orbit case 



The relevance of the circular orbit case is supported by the fact that gravitational radiation can circularize 
the motion and drive the binary toward the innermost stable circular orbit. Moreover, the dynamics can be 
integrated explicitly in time in this limit. Although the main steps of the method given above do not change 
some equations and expressions become simpler. 

The relative velocity vector can be generally decomposed as 



V = rn + rcjm 



(24) 



and circular orbits are defined by the f = and lu — conditions. To obtain the expressions for the polarization 
states we use the relevant terms in the decomposition of the relative velocity given in Ref. [E] : 



2 2 22 

V = V I = r UJ = 



M 



l-]f^(2 + 3C)(nxm)S 



3/2' 



Then the contributions for and ft,x which change compared to the general formula become 

M 



l,bSO 



= [- [-{Pi-ql) + ipt~ql)] 



2v 



[Sz{pI ~ ql) + 2Sx{pxPz - qxqz)] 



+ 72 ['^^(pI " 9'^) + '^SzipxPy - qxqy) + 5Sx{PxPz - qxqz)+ 
+ 2Sx{pyPz - qyqz)-2Ny{{qS)p^ + {pS)q^)-2N^{iqS)py + {pS)qy)] 



and 



1,550 



(~f) i~P^1^ +Pvqy] > 



n- 



+ 



vC 



[SzPxqx + SxiPxqz + qxPz)] + 

[2SzPxqx - 2Sz{pxqy + qxPy) + ^Sx{Pxqz + qxPz 



'iSxiPyqz + qyPz)~ 2Ny{{ciS)q^ ~ {pS)p^) -2N^{{ciS)qy - {pS)py 



(25) 



(26) 



(27) 
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5 Extreme mass ratio limit 



An important special case discussed in the literature is the extreme mass ratio limit, when the mass ratio of 
the bodies is negligible. We consider the case when the mass of the rotating body is much greater than the 
other. This way 77 « ^ <C 1, A/ « mi, /i w m2 and 5m/m k, 1 and the description of the motion and the formal 
expressions for the polarization states change. 

In the extreme mass ratio limit the components of the relative velocity become 



2E 2mi 

7712 r 

L 



2E 

TO2 



2toi 



2„2 



4LS 

m2r^ 



v± 



Lrs 
m2r% 



2LS 



m2rN 

and the equations of the motion for the angle variables are: 

L 



is sin $7v 
Ls cos $ N 



2J 



— T 



N 



' N 



m2r^ 



2Lrs 2LS 



1712 



cosS, 



m2 



cos /Sat — sill PnPpn 



sinaS" 



"^2 . > m2 

sm^ sm Pat + cos PnPpn 



L 



The contributions for the polarization states can be written as 



hi = 



"^1 w . 



:) + 2vj_r{p^py ~ q^qy) + v^ipl - ql) 



l,bSO 



4" {iVA.Szipl 



(pS)<z. 



(28) 



(29) 



-r[S ^{pxV-qx(\)\x-2v±_[S X (Pa:P-9xq)]y} , 



(30) 



and similarly 



ft^ = 2 



l,5SO 



^2 _ TOi 



Pxq, 



+ v±r{pa;qy + qxPy) + vlPyqyj , 



= ^ [(qS)gx - {pS)px] , 
2 

= {Qv±SzPxqx + r[S x{p^(i+qxp)]x-2v±[S x {pxCi+qxp)]y} . 



(31) 



In a previous work |17| the polarization states was described in the Lense-Thirring approximation. The z 
axis of the invariant system was fixed to S since the precession of the spin is negligible in that case. There 
is a constant rotational transformation between the above description and the Lense-Thirring case. The main 
difference is that in the Lense-Thirring approximation J fails to be zero. If one specifies that J = in this 
comparable mass case the common limit of the two descriptions can be found. 



6 Conclusions and remarks 

In this article we have presented a method to evaluate the detectable gravitational wave signals generated by 
a spinning compact binary system moving on eccentric orbit in the case of one spinning object up to 1.5 PN 
order. We have introduced the invariant and a comoving coordinate systems to describe the evolution of the 
dynamical quantities. With the use of the constants of motion we have discussed the equations describing the 
evolution of the dynamical quantities for the determination of the polarization states. We have calculated the 
components of the relative velocity vector, the spin and the (N,p,q) triad in terms of these quantities, namely 
the length of the separation vector and the Euler-angles. To determine the effect of the eccentricity of the orbit 
on the detectable signals we have investigated the circular orbit limit. An other significant property of this case 
is that the explicit time dependence of /i+ and hx can be calculated. 

The results presented here are independent of the parameterization of the orbit. We plan to use the general- 
ized true anomaly parameterization of the motion [19 to investigate the structure of the wave signals. Moreover, 
this method can be the starting point to study the properties of two spinning objects and unbound orbits [22j . 
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A The transverse-traceless tensor 

The decomposition of the relative velocity vector, Eq. results the following form of the quadrupole and 
spin-orbit terms of hl^j^: 



Q'^ = 2 



, M 




PQlo = --^[SxN](V), (33) 



^'r^•'' + 2v±rm'^'n^'> + v^mj^m^ , (32) 

;■ J 

2M 

P^'^Q'^^ = ( 6v±M\n X rnlSn'n^ + (2mi - m2)fn^'\n x Sl^^ 

— (5m2 + 4mi)w^n^*[m x S]-'-' — 2Tn2V±m^^[m x S]-'^ 

- {fN^ + 2v±Ny)m2[S X Nfn^^ - 2N^v±m2[S x N](*m^')) . (34) 
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